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Abstract: We perform a comprehensive study of the impact of new-physics op- 
erators with different Lorentz structures on B — > D* + l~vn decays, t = e~,fi~,r~, 
involving the b — > clvg transition. We present the full three angle and q 2 angular 
distribution with new physics operators with complex couplings. Various observ- 
ables are constructed from the angular distribution and we specially focus on the CP 
violating triple product asymmetries which vanish in the Standard Model without 
any hadronic complications. Two of the three triple products are only sensitive to 
vector/axial vector new physics and do not depend on pseudoscalar new physics. 
Hence the measurements of non zero triple product asymmetries will be clear sign of 
new physics and strong signals for vector/axial vector new physics operators. Even 
though we focus on r final state, one can use the triple products to search for new 
physics with e and [i in the final state. 
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1. Introduction 

The Standard Model (SM) of particle physics, even though very successful, is ex- 
pected to break down at some energy scale and make way for a more complete 
theory. Exploration of what lies beyond the SM can be carried out at the energy 
frontier in colliders such as the LHC or at the intensity frontier at high luminosity 
experiments. In the intensity frontier, the B factories, BABAR and Belle, have pro- 
duced an enormous quantity of data in the last decade. There is still a lot of data 
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to be analyzed from both experiments. The B factories have firmly established the 
CKM mechanism as the leading order contributor to CP violating phenomena in the 
flavor sector involving quarks. New physics (NP) effects can add to the leading order 
term producing deviations from the SM predictions. In this respect, the second and 
third generation quarks and leptons are quite special because they are comparatively 
heavier and are expected to be relatively more sensitive to NP. As an example, in 
certain versions of the two Higgs doublet models (2HDM) the couplings of the new 
Higgs bosons are proportional to the masses and so NP effects are more pronounced 
for the heavier generations. Moreover, the constraints on NP involving, specially the 
third generation leptons and quarks, are somewhat weaker allowing for larger NP 
effects. 



It is interesting that there are certain discrepancies in decays involving r and v T 
states though none of them are significant enough to establish clearly the presence 
of NP. There is a seeming violation of universality in the tau lepton coupling to the 
W suggested by the Lep II data which could indicate NP associated with the third 
generation lepton 0. Recent measurement of CP violation in r decay find Aqp 
in t~ ir~K s {> Otv°)u t is (-0.36 ± 0.23 ± 0.11) % which is different from the SM 
prediction (0.36 ± 0.01) by 2.8 a. The branching ratio of B — > tv t showed some 
tension with the SM predictions || and this could indicate NP, possibly coming 
from an extended scalar or gauge sector M. However, new Belle || and BaBar 
II measurements, obtained using the hadronic tagging method, are more consistent 
with the SM. 



If there is NP involving the third generation leptons one can search for it in 
semileptonic b —> ctu t decays such as B — >■ D + t~i7 t , B — > D* + t~v t [[/J. The 
semileptonic decays of B meson to the r lepton is mediated by a W boson in the SM 
and it is quite well understood theoretically. In many models of NP this decay gets 
contributions from additional states like new vector bosons or new scalar particles. 
The exclusive decays B — > D + t~i7 t and B — > D* + t~i7 t are important places to 
look for NP because, being three body decays, they offer a host of observables in 
the angular distributions of the final state particles. The theoretical uncertainties 
of the SM predictions have gone down significantly in recent years because of the 
developments in heavy-quark effective theory (HQET). The experimental situation 
has also improved a lot since the first observation of the decay B — > D* + t~i7 t in 2007 
by the Belle Collaboration ||. After 2007 many improved measurements have been 
reported by both the BaBaR and Belle collaborations and the evidence for the decay 
B — > D + t~i7 t has also been found || [10], [lTJ. Recently, the BaBar collaboration 
with their full data sample of an integrated luminosity 426 fb _1 has reported the 
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measurements of the quantities ||T2 

BR(B D + t-j7 t ) 



R{D) 
R(D*) 



BR(B -> D+£-u e ) 
BR(B -»■ J D*+r-z7 r ) 



The SM predictions for and R(D*) are 0, 0, [14 

R(D) = 0.297 ±0.017, 
R(D*) = 0.252 ±0.003, 



0.440 ± 0.058 ± 0.042 , 
0.332 ±0.024 ±0.018 



;i-2) 



which deviate from the BaBar measurements by 2cr and 2.7a respectively. The 
BaBar collaboration themselves reported a 3. 4a deviation from SM when the two 
measurements of Eq. (|1.1| ) are taken together. 

These deviations could be sign of NP and already certain models of NP have 
been considered to explain the data O, IT5L [1 



mm 



gg, grj, 0, |23], g, fzg. 

In Ref. fTTfl , we calculated various observables in i? — y D + t~i7 t and i? — y D* + t~i7 t 
decays with NP using an effective Lagrangian approach. The Lagrangian contains 
two quarks and two leptons scalar, pseudoscalar, vector, axial vector and tensor 
operators. Considering subsets of the NP operators at a time, the coefficient of 
these operators can be fixed from the BaBar measurements and then one can study 
the effect of these operators on the various observables. In this work we extend 
the work of Ref. JT7J by providing the full angular distribution with NP. The full 
angular distribution, in the SM, has already been used in experiments for final states 
with n and the electron. We are proposing to study the angular distribution with 
r final states. We will also discuss CP violating observables which are the triple 
product (TP) asymmetries [27]. In the SM these TPs rigorously vanish and so any 
non-zero measurements of these terms are clear signs of NP without any hadronic 
uncertainties. Moreover, as we will see most of the TPs depend on the vector/axial 
vector couplings and not on the pseudoscalar couplings. Hence these TPs provide 
useful clues to the nature of NP. As in the previous work we will neglect the tensor 
term in the effective Lagrangian. 

The paper is organized in the following manner. In Sec. 2 we set up our formalism 
where we introduce the effective Lagrangian for NP and define the various observables 
in B — y D* + t~i7 t decays. We also present the SM predictions for these observables in 
that section. In Sec. 3 we present the numerical predictions which include constraints 
on the NP couplings as well as predictions for the various observables with NP. 
Finally, in Sec. 4 we summarize the results of our analysis. 

2. Formalism 



In the presence of NP, the effective Hamiltonian for the quark-level transition b — y 
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cl v\ can be written in the form ]28 
AG F V cb 



H, 



eff 



V2 



1 + V L ) [^P L b] [l-fP LVl ] + V R [crfP R b] [h»P L n] 



+S L [cP L b][lP L v l ] + S R [cP R b] [lP m ] +T L [ca^P L b] [Ta^Pm] (2.1) 



where Gf = 1.116637 x 10~ 5 GeV~ 2 is the Fermi coupling constant, V c b is the 
Cabibbo-Koboyashi-Maskawa (CKM) matrix element. Pl.r — (1 T 7s)/2 is the 
projector of negative/positive chirality, and we use cr^ = i[ 7At , 7i/]/2. We have as- 
sumed the neutrinos to be always left chiral. Further, we do not assume any relation 
between b — > ul~vi and b — > cPv\ transitions and hence do not include constraints 
from B — > tv t . The SM effective Hamiltonian corresponds to Vl — Vr — Sl — 
S R = Tl = T R = 0. In this paper we will ignore the tensor interactions. With this 
simplification we write the effective Lagrangian as 



H, 



eff 



GpVcb 

~7T 



C7^(l - 75)6 + gycy^b + g A cj^ 5 b 1^(1 ~ 7s)^ 



+ 



g s cb + g P cry s b 1(1 - 75)^ + h.c 



(2.2) 



where g ViA = V R ±V L and g s ,p = S R ±S L . 

The effects of NP can be seen in the helicity amplitudes that describe the decays. 
The expressions for the hadronic helicity amplitudes for the B — > D*tu t decays are 



2rriD* v q 2 



A 



(m| - m 2 D * - q 2 )(m B + m D .)Ax(g 2 ) 

A\\ = V2(m B + m D *)^i(g 2 )(i - gA) , 
A, 



Am 2 B \p D *\ 2 
m B + m D * 



A 2 (q 2 ) (l-g A ) 



-V2- : -\p D *\{l + g v ) 



(m B + m D * 
^ = 2m B \p D *\A (q 2 ) 



9 a) 



A 



2m B \p D *\A (q 2 



:9p 



(m b (fi) + m c (/i))' 
where the t and the P amplitudes arise in the combination 



(2.3) 



A, 



tp 



A t 



-A, 



(2.4) 



The form factors Ai 2 fi(q 2 ) and V(q 2 ) are defined in the appendix. As is clear from 
the above equation various helicity amplitudes are sensitive to different NP operators. 
These helicity amplitudes can be probed in various differential distributions providing 
useful information about NP. 
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The transversity amplitudes A\\ and A± are related to the helicity amplitudes 
A± as 

A± = 71 (Af ~ AJ) ' 

4, = (.A+ + .4.) . (2.5) 

All the amplitudes are complex if the NP couplings are complex. The phases in 
the couplings are weak phases and change sign when we go from particle to anti- 
particle decays. Though strong phases in the current can arise from higher order 
loops these will be tiny and we will ignore them. Hence the only CP violating signals 
will be of the triple product types and all direct CP violating effects will vanish. 
Moreover we see that Ao t \\,t have the same weak phases and any interference between 
these amplitudes will not lead to any CP violating signals. The only CP violating 
signals will come from the interference of A± with the other vector/axial vector and 
pseudoscalar current amplitudes ^4o,||,tP- 
We will now consider the two cases: 

• Case a : In this case we will set Sl, Sr = and assume that the NP affects lep- 
tons of only the third generation. This scenario could arise from the exchange 



of a new charged W boson p9fl . We point out that this is just a simplifying 
assumption and in fact the general angular distribution presented in the paper 
is applicable to also e and \i in the final state. 

Case b : In this case we will set Vl, Vr = and assume that the NP affects 
only leptons of the third generation. This scenario could arise in models with 



extended scalar sector [30 



2.1 Constraints on NP couplings 

For the numerical calculation, we use the B — > D and B — > D* form factors in the 
heavy quark effective theory framework |31], |32| . B — > D* form factors are summa- 
rized in appendix. The constraints on the complex NP couplings in the b — > cl~V[ 
effective Hamiltonian come from the the measured R(D) and R(D*) in Eq. ( |1 . 1| ) at 
95% C.L. We also vary the free parameters in the form factors discussed in appendix 
within their error bars. All the other numerical values are taken from [RHl and 



A detail analysis of R{D) and NP in the decay B — > D + £~i>£ can be found in [17]. 
The allowed ranges for NP couplings are then used for predicting the allowed ranges 
for the observables in the the angular distribution discussed below. 

The combination of the couplings gv = Vr + Vl appears in both R(D) and 
R(D*), while Qa — Vr — Vl appears only in R(D*). Vr and Vl receive constraints 
from both R(D) and R(D*). If NP is established in both R(D) and R(D*) then the 
case of pure coupling is ruled out. The constraints on the complex couplings gy 
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and qa are shown in the colored region of Fig. [I] (left) and (right). We confirm from 
Eq. Q2.3|) that if the new interaction in purely left-handed then the amplitude and all 
the distributions just gets scaled by a common factor. Hence instead of considering 
the pure V — A and V + A quark current cases we will consider cases which include 
pure gv or pure qa complex couplings. Interestingly, the analysis in Ref. |HJ indicates 
that the data prefers either pure vector or pure axial vector couplings. 



Only g v presents 
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Figure 1: The constraints on the complex coupling gy = Vr + Vl (left panel) and gA = 
Vr - Vl (right panel) at 95% C.L. 

The combination of the couplings g$ — Sr + Sl appears only in R(D), while 
gp = Sr — Sl appears only in R(D*). If NP is established in both R(D) and R(D*) 
then the cases of pure gs or gp couplings will be ruled out. The constraints on the 
complex couplings Sr + Sl and Sr — Sl are shown in Fig. ||. The constraint on the 
complex Sr and Sl from the recent R{D) and R(D*) measurements are relatively 
weaker than those in the (Vl, Vr) case. 
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Figure 2: The constraints on the couplings for complex (Sr + Sl) (left panel), complex 
(Sr - S L ) (right panel) at 95% C.L. . 
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3. Angular analysis 

The complete three-angle distribution for the decay B — > D*{— > Dtt)1^i?i in the 
presence of NP can be expressed in terms of four kinematic variables q 2 , two polar 
angles 9i, 9d*, and the azimuthal angle x- The angle 9 t is the polar angle between the 
charged lepton and the direction opposite to the D* meson in the (Ivi) rest frame. 
The angle 9d* is the polar angle between the D meson and the direction D* meson 
in the {Dir) rest frame. The angle x is t ne azimuthal angle between the two decay 
planes spanned by the 3-momenta of the (Dir) and (Ivi) systems. These angles are 
described in Fig. 0. The three-angle distribution can be obtained by using the helicity 
formalism. 

We can write the angular distribution explicitly for easy comparison with previous 
literature fl35l B6l B7I, M 



9 



dq 2 d cos 9[ d cos 9 D * dx 327r 



i=l H j=l 



(3.1) 



where 



h 
Ji 

h 

■h 
h 
■h 
h 

h 

h 
■h 
h 
■h 
h 



4 cos 2 9 D * sin 2 9i\A \ 



4 cos 2 9 D , 



|A| 2 cos 2 6»i + \A tP \ 2 - 2Re[A tP A* ] cos 6^ 



sin 2 # 



D* 



{\A\\\ 2 + \A±_\ 2 ){1 + cos 2 9i) - 4Re[A\\A* ± ] cos9 t 



sin 2 9 D * sin 2 9i{\A\\\ 2 + \A_ 



I 2 ) 

|2 



|Al| 
|2\ 



— sin 2 9d* sin 2 Q\ cos 2x(|.A|| 
sin 2 ^* sin 2 9i cos 2x(\A\\ \ 2 - \ A±\ 
-2 v / 2sin2^ D * sin 9i cos xRe[A±Al] , 

sin 29 d* sin 9 h cos xRe[A\\A* tP \ , 

sin 9\ cos 9\ cos xRe\A^A*§ 
— 2\l~2 sin 29b* sin 9i cos9i cos xRe[A\\A%\, 
sin 2 9d* sin 2 9\ sin , 

— sin 2 6 D * sin 2 61 sin 2xIm[A\\A 



l 0J) 



-2v / 2sin2^ D . sin 9i sin x/m[y4.||.4j||], 
-2V2 sin 29 d* sin 9i sin xIm[A_LA* P ] , 

sin 29\ sin % Jm [.Aj^q] , 
- V2 sin 29 D * sin 29i sin [Al-Aq] > 



where the quantity N F is 

-G\\v 



D* 



IT/ \ 2 n 2 

\ v cb\ Q 



3 x 2 6 vr 3 m| 



Til ^ \ 2 

Br(D*^D7i) 

q 2 J 



(3.2) 



(3.3) 
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Figure 3: The description of the angles Q\ y D* and X i n the angular distribution of B 
£)*(->■ Dn)l~vi decay. 



When there are no strong phases then A\\ and Ao have the same weak phase and If 
vanishes. The momentum of the D* meson in the B meson rest frame is denoted as 
\Vd*\ = A 1 / 2 (m^, m 2 D *,q 2 )/2m B with A(a, b, c) = a 2 + b 2 + c 2 — 2(ab + be + ca). 

The complex couplings leads to CP violation which are sensitive to the angular 
terms sin x and sin 2%. The coefficients of these terms are TPs and have the structure 



Im[AiA*} ~ sin( 



bj), where Aij 



\Ai t j 1 6 



i<t>i. 



In the SM these terms vanish 



as there is only one contribution to the decay and so all amplitudes have the same 
weak phase. For the charged conjugate modes the weak phases change sign and 



A 



M 



\A 



Ml 



-i<t>i, 



and the TPs change sign. Hence any non-zero measurements of 
the TPs is a clear sign of NP without any hadronic uncertainties. Even though we 
focus on t final states, we should point out that this distribution is applicable also for 
e and fi in the final state. Since experiments have already studied this distribution 
for e, fj, final states it might be worth checking the sinx and sin 2% terms in the 
distributions for these decays for signals of non SM physics. 

It will be convenient to rewrite the angular distribution as p9| , 



d 4 T 



dq 2 d cos 9i d cos 9 b* d\ 



9 

— NF 
32vr 



cos 2 9 D * {v? + 1/ 2 ° cos 26 1 + V£ cos 9 t ) 

+ sin 2 9 D . (V? + V? cos 20, + V? cos 0j) 

+ \/ 4 T sin 2 9 D * sin 6i cos sin 26i cos x 

sin 9 1 cos x + sin 2 9d* sin 2 9i sin 2\ 

+ V° T sin 29 D * sin 9 t sin X + V° T sin 29 D « sin 29] sir 



The decay B — > £)*(—>• Dtt)1 v\ is completely described in terms of twelve angular co- 
efficient functions V*. These angular coefficients depend on the couplings, kinematic 
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variables and form factors, and are given in the Appendix in Eq. ( B.7 ), Eq. (| 
and Eq. ( |B.9|) . The definitions of the angles involve the directions of the l~ and D*. 

For the CP-conjugate decay B — > £)*(—>• Dir)l + vi, one defines these angles rela- 
tive to the directions of the r + and D*. The V^'s can be obtained from the V^'s by 
replacing 9\ —¥ 9i + 7r and x ~Xi an d changing the signs of the weak phases. This 



transformation is equivalent to replacing V® 2 ~~ >" Vi2> ^3° ~~ ~ ^3°; ^i2i ~^ V\ 



2,4' 



— fj = ^Nf\ cos 2 9 d « (V? + V 2 ° cos 20, - V 3 ° cos d l 

dq d cos 9\ d cos 9tj* d\ o2ir I V 

+ sin 2 9 D , (W + V 2 T cos 29 l - V£ cos { ) 

+ sin 2 9 D * sin 2 9 X cos 2\ + V° T sin 20 D . sin 2Q X cos x 
— V 2 T sin 26 , £ ) * sin 0/ cos x — sin 2 9 D * sin 2 0; sin 2% 

+ V° T sin 20£>* sin 0, sin x - sin 20£>« sin 20; sin x 

(3.5) 

3.1 Differential branching ratio 

The angular distribution allows us to define several observables. The starting point 
is to obtain the decay rates dT /dq 2 after performing integration over all the angles 

dT 8irN F . A 
- = — { A L + A T ), 

(3.6) 

where the D* longitudinal and transverse polarization amplitudes Al and At are 
A l=( v ?-\v?)> A T = 2(V?- 1 ~V?). (3.7) 



One can see from Eq. ( |B.7|) and Eq. (p.8j) that A L is proportional to \Aq\ and |^4<p| 
while At to |-4.||| 2 + |-4_l| 2 - The D* polarization amplitudes can be extracted from 
the angular distribution in cos9d* (see Eq. (|3.10|) below). Since there are no direct 
CP violation we have Al t — Al,t- Hence, 

dT _ dT 

dq 2 ~ dq 2 ■ {6 * ] 

Furthermore, one can also explore the q 2 dependent ratio 

dBr[B -» D*+T-v T }/dq 2 
W ) = d Br[B ^ DH-^/dq 2 ' (3 ' 9) 
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9 2 [Gev 2 ] 9 2 [G«v 2 ] 

Figure 4: The left (right) panels of the figure show DBR (Rp,*{q 2 )) for the decay B° — > 
D* + tv t in the scenario where only g A (upper) and only gp (lower) couplings are present. 
The green band corresponds to the SM prediction and its uncertainties; The yellow and 
purple bands show predictions for allowed values of only gA and only gp, respectively 

where I denotes the light lepton (e,/z). The ratio Rp>* are independent of the form 
factor h Al {w). 

Fig. [| shows DBR and Rp,*(q 2 ) for B° — > D* + tv t in the presence of only g A = 
Vr — Vl and only gp = Sp — Sl- We make the following observations: 

• If only g A couplings are present, DBR can be enhanced up to 0.4% at q 2 ~ 
8.5GeV 2 . Ro*{q 2 ) can be enhanced up to 0.9% at high q 2 . 

• If only gp couplings are present, DBR can be enhanced up to 0.4% at q 2 ~ 
7.5GeV 2 . Note that the peak of the DBR is shifted to low q 2 direction relative 
to the SM. Rp,*(q 2 ) can take the value of about 0.7 at q 2 m 7.5GeV 2 . 

Finally, the new NP coupling gy only appears in the transverse amplitude A±, 
and does not significantly affect the DBR and Rp,*(q 2 ). 

3.2 Polarization fraction for D* 

The differential angular distribution in cos dp,* gives access to the polarization frac- 
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tion of the D* meson in the decay B — > £)*(—>■ Dn] 



T U-r 



d 2 T 
dq 2 d cos 9d> 



1 dV 



D * cos 2 



e D * + (i-Fr. 



sin 2 9 D * 



(3.10) 



where we define the longitudinal and transverse polarization fractions of the D* 
meson as 



A L + A T 



A L + A T 



(3.11) 



with Ff 



+ FP* 



T — 1. Similarly, one can define the polarization fractions Fp* T for CP- 
conjugate mode but they are the same as Fp T in the absence of direct CP violation. 

When only the coupling is present the polarization fractions of the D* meson 
gets contributions from the amplitudes Ao and A\\ , which are functions of the new 
coupling. Due to cancellation of the Qa coupling contributions in Fj? , its behavior 
is similar to its SM prediction. New gv coupling appears only in the amplitude A± 
and again F®* looks similar to its SM prediction. In the plot Fig. |5] we show F£* in 
the presence of the gp coupling. In this case F£* can be as large as 0.85 at low q 2 , 
and it decreases to the SM value at high q 2 . 
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Figure 5: The figure shows Fp* for the decay B — > D* + t~i7 t in the scenario where 
only gp coupling is present. The green band corresponds to the SM prediction and its 
uncertainties; The purple band shows prediction for allowed values of only gp. 



3.3 Distribution in cos#/ and App 

The forward-backward asymmetry App can be obtained from the single-differential 
angular distribution: 



d 2 Y 
dq 2 d cos 9 1 



87rN f 



(Vi T + \v?) + (Vf + \vf) cos 20, + (Vf + l -V°) cosOi 



(3.12) 
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The forward-backward asymmetry (FBA) for the leptons is defined by 

L d cos 81 , 2 d j T a — f -, d cos 9i -i 2 f r a 

a I 2\ J U ' dq A a cos J— 1 l dq z a cos Oi /q i q\ 

Jo d C0S g < ^toflj + J-l d CQS ^ d^dcose, 

Then one can obtain: 

similarly, FBA for the conjugate mode is 

r V T 4- 

>w)=- ■ (3 - i5) 

If the absence of direct CP violation A F B(q 2 ) = —A F B.(q 2 )- We define the average 
FBA as : 

lA FB } D *{q 2 ) = \{A FB {q 2 ) - A FB {q 2 )) . (3.16) 

Within the SM, [A FB ] D * (q 2 ) has a zero crossing at q 2 « 5.64GeV 2 (see Fig. (|])). 
We make the following observations from this figure: 

• If only the Qa or only the gy couplings are present, the FBA can reach a value 
close to 50% at low q 2 and its sign is mostly negative. The FBA converge to 
its SM prediction at high q 2 . 

• If only the gp coupling is present, the FBA can reach a value up to 30% at low 
q 2 . It can have both positive or negative signs. Again, the FBA converge to 
its SM prediction at high q 2 . 

3.4 Asymmetries Aq in the angular variable x 

In this section, we consider three different transverse asymmetry A^ -1 ' 2 . These 
asymmetries are obtained by integrating out the polar angles #/ and Q B * in the three 
different regions. 

3.4.1 A$ 

The transverse asymmetry A^j is defined through the angular distribution in x as 



d 2 T 1 dT 



1 + (Aj. 1) cos2x + 4! ) sin2x) . (3.17) 



dq 2 d\ 27r dq 2 

It can be obtained by integrating Eq. (|3.4j) over the two polar angles d\ and Or,*. 
Here A^ is a TP, and is discussed separately below. One can obtain 
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Figure 6: The figures show L4_fb]d* for the decay B — > D* + t~i7 t in the scenario where 
only gj\ , only gy and only gp couplings are present. The green band corresponds to the 
SM prediction and its uncertainties; The yellow, red and purple bands show predictions 
for allowed values of only gA, only gv and only gp, respectively. 



similarly, for the conjugate mode is 



4V) = ,t 4 t (3-19) 
In the absence of direct CP violation Aq = Aq . We define the average A^(q 2 ) as 

<4V)) = ^(4V)+4V))- (3.20) 

The SM prediction for (A^\q 2 )) is shown in the green band of Fig. [7[ Similar to 
FjP* , the asymmetry (4P) remains almost the same as the SM prediction when only 
gA or only g v couplings are presents. 

New gp coupling appears only in the amplitude Ap and affects only the denom- 
inator of Aq\ The purple band in Fig. [7] is the prediction for (4? ) when only the 
gp coupling present. It cannot be enhanced more than 5% in this scenario. The 
amplitude Ap vanishes at high g 2 , and hence (4? ) reduces to its SM value. 



3.4.2 A c 



(2) 
C 

We define the angular distribution 
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Figure 7: The figure shows (A c '(q 2 )) for the decay B — > D* + t~i7 t in the scenario where 
only the gp coupling is present. The green band corresponds to the SM prediction and its 
uncertainties; The purple band shows prediction for allowed values of only gp. 



d 2 rw 

dq 2 dx 
One can obtain: 



o ., 



d 4 r 



_i J dq 2 d cos 6i d cos Op,* d\ 



d cos9p* 



(3.21) 



d 2 T^ 1 dT 



dq 2 dx 4 dq 2 



A c ' cosx + A T ' sm X 



(3.22) 



where 



(2)^2> 



(A L + A T ) 



(3.23) 



(2) 

Here A T ' is a TP, and is discussed separately below. For the conjugate mode 



(2)^2> 



-vr 



(A L + A T ) 



(3.24) 



In the absence of direct CP violation A 



$ = —Aq. We define the average A c V) (q 2 ) 



as 



(4V)> = ^(W)-W))- 



l(2)/ 2^ 



(3.25) 



(/^(g 2 )) depend on all three new couplings gy and gp. For all g 2 , the magnitude 
of (Aq) is generally suppressed by these new couplings. As shown in Fig. [8], in all 
three cases, the value of (Aq ) can be anywhere in the q 2 range. It can be either 
positive or negative. In particular, there may or may not be a non-SM zero crossing. 
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Figure 8: The figures show (A^\q 2 )} for the decay B° — > D* + tu t in the scenario where 
only gA , only gy and only gp couplings are present. The green band corresponds to the 
SM prediction and its uncertainties; The yellow, red and purple bands show predictions 
for allowed values of only gA, only gy and only gp, respectively. 



3.4.3 A 



(3) 
C 



Finally, we define the single angle distribution 



d 2 T^ 
dq 2 dx 



1 r 



J-l 



d cos 9i 



1 r . 



J-l 



d cos Op,* 



d 4 r 



dq 2 d cos 6 1 d cos 9p* d\ 

(3.26) 



One can obtain: 



d 2 r^ _ 2_ dr_ 

dq 2 dx 37r dq 2 



a (3) . .(3) . 

A c cosx + Arp smx 



(3.27) 



where 



or 



(A L + A T ) 

(3) 

Here is a TP, and is discussed separately below. For the conjugate mode 

Vi 0T 



(3.28) 



(3.29) 
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In the absence of direct CP violation A c ' = A c . We define the average A c (q 2 ) as 

(4 3) (g 2 )) = ^(4 3) (^) + 4 3) (^))- (3-30) 

The angular coefficients V® T depends only on the coupling. Due to cancella- 
tions of the NP contributions, (A^q) behaves similar to its SM prediction when only 
the gA coupling is present. The SM prediction of (A^\q 2 )) is shown in the green 
band of Fig. [|. (Aq^) is not sensitive to the new <jy coupling. As shown in the purple 
band in Fig. |9], (Aq) is suppressed relative to the SM by the new gp coupling. 

Only g P presents 

0.20 
0.15 

A 

~y 0.10 

V 

0.05 
0.00 

4 5 6 7 8 9 10 
g 2 [GeV 2 ] 

Figure 9: The figure shows (AP(q 2 )) for the decay B° -> D*+ tv t in the scenario where 
only the gp coupling is present. The green band corresponds to the SM prediction and its 
uncertainties; The purple band shows prediction for allowed values of only gp. 




3.5 CP- violating triple-product asymmetries 

In this subsection, we consider the TPs in the decays B — > D*{— > Dii)l~&i and B — > 
D*(—> Dir)l + i>i. For the decaying B meson, the TP is proportional to {hp, x hi) ■ h z 
in its rest frame, where the unit vectors are given in terms of the momenta of the 
final-state particles as [[39 



n D = r-~ —\i n z = r« — i ' r = i > ' i ' n i = r« • ( 3 -3i) 

\Pd x p w \ \Pd+pA \Pi- x pp T \ 

The vectors hp and hi are perpendicular to the decay planes of the D* and the 
virtual vector boson. In terms of the azimuthal angle \i one gets 

cosx = h D -hi, smx = (h D x hi) ■ h z , (3.32) 

and hence the quantities that are coefficients of sin x (or of sin 2\ = 2 sin \ cos x) are 
the TPs. 

As noted above, while the angular distribution for the B decay involves x, for B 
it involves — X- 
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Figure 10: The left (right) panels of the figure show (Alp) for the decay B — > D* + t v t 
in the scenario where only qa (yellow band) and only gy (red band) couplings are present. 



3.5.1 A 



U) 



The first TP is Alp, introduced above in Eq. (|3.17|) . One can find A ( r> and A ( p as 



(i) 



4V) 



W * A?( q 2 )-- 



3(A L + A T 



3(A L + A T ) 



(3.33) 



In the absence of direct CP violation AP = Asp. We observe that Ap depends 
on both qa and gy couplings and not on the gp coupling. The CP- violating triple- 
product asymmetry is 



<4V)> = ^(Aipur) , - APuh 



(3.34) 



Fig. |10| shows (A$p) for B — > D* + t v t in the presence of only qa and only gy 
couplings. We make the following observations: 

• If only the qa coupling is present, the magnitude of (Ajp) can be enhanced up 
to 4% at q 2 ~ 8.0GeV 2 . It vanishes at the end points as the amplitude A± 
diminishes. (AP) can have either positive or negative. It may or may not has 
non-SM zero crossing points. 

• If only the gy couplings is present, the magnitude of (AP) can be enhanced 
up to 5% at q 2 m 8.0GeV 2 . The behavior of (Aip) is similar to the above case. 



3.5.2 A 



(2) 
T 



The second TP is A$p , introduced above in Eq. ( [3.221 ). AP and A$p are given by 

4V) = ttTtv 4 2) = 7T^rr- (3-35) 



(A L + A T ) 



(A L + A T ) 
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We observe that A T depends on all three new couplings g& , gy and gp. This TP is 
proportional to the lepton mass and so is very small when the lepton is the electron 
or the muon. The CP-violating triple-product asymmetry is 

<4V)) = ~(4V)-4V))- (3-36) 



Fig. |TT] shows (A T ) for B — > D* + t~v t in the presence of only g& , only g v and only 



gp couplings. We make the following observations: 

• If only the g& coupling is present, the magnitude of (AP) can go up to 10% at 
low q 2 . (A$p) vanishes in the end points as the amplitudes A± and A t p fall. 
It can have either sign at both low and high q 2 . Also (4^) may or may not 
has non-SM zero crossing. 

• If only the gy coupling is present, the magnitude of (AP) can reach up to 
10% at low q 2 . The behavior of (AP) is similar to the one when only the g^ 
coupling is present. 

• If only the gp coupling is present, the asymmetry prediction is similar to the 
other two cases. 



3.5.3 AP 



The third TP is introduced above in Eq. (|3.27 ). Ap and AP are given by 

4V)= ,/ 4 A r , 4 3) = ~ ,/ 4 t , • (3-37) 
(A L + A T ) t (A l + A t ) 

(3) 

We observe that A T depends on both the new couplings g& and gv but does not 
depend on g P . The CP- violating triple-product asymmetry is 

(A?( q 2 )) = l(AP(q 2 ) + AP(q 2 )). (3.38) 

Fig. 12 shows (AP} for B — > D* + t~i7 t in the presence of only g^ and only gy 
couplings. We make the following observations: 

• If only the coupling is present, the magnitude of {Ap) can be enhanced up 
to 4% at q 2 ~ 8.0GeV 2 . It vanishes at the end points as the amplitude A± 
diminishes. (AP) can have either sign at both low and high q 2 . Also it may 
or may not has non-SM zero crossing. 

• If only the gy coupling is present, the magnitude of (AP) can be enhanced up 
to 5% at q 2 fa 8.0GeV 2 . The behavior of [Aj, ) is similar to the case above. 
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Figure 11: The left, center and right panels of the figure show (Aj, ) for the decay 
B — > D* + t~i7 t in the scenario where only gA (yellow band), only gy (red band), and only 
gp (purple band) couplings are present. 
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Figure 12: The left (right) panels of the figure show (^4^ ) for the decay B — > D* + t~i7 t 
in the scenario where only gA (yellow band) and only gv (red band) couplings are present. 



4. Discussion and Summary 

In summary we presented the full three angle and q 2 distribution for B — > D*l~vg 
though our focus was on the decay B — > D* + t~i7 t because of new experimental results 
that are not consistent with SM predictions. We extended the work of Ref. [17| by 
considering additional observables from the angular distribution. Particular attention 
was paid to the CP violating triple product asymmetries. It was argued that in the 
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SM these asymmetries vanish rigorously and so a non zero measurement would be 
a smoking gun signal for new physics. Of the three triple product asymmetries two 
are sensitive to only vector and axial vector new physics. Hence the triple product 
asymmetries are not only sensitive to new physics but also can probe the nature of 
new physics. 
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A. Details of the B — > D*r u r angular analysis 
A.l kinematics 

In the B rest frame, the co-ordinates are chosen such that the D* meson is moving 
along the positive z-axis, whereas the virtual gauge boson is moving along the nega- 
tive z-axis. The four-momenta of the B and D* mesons, and the virtual gauge boson 
are 

Pb = (m B , 0, 0, 0) , p D * = (E D *,0,0, \p D *\) , q = (q , 0, 0, -\p D *\) , (A.l) 

where E D * = {m\ + m 2 D , — q 2 )/2m B , \pd*\ = A 1//2 (m^, m 2 D *, q 2 )/2m B , and q = 
(m 2 B — m 2 D * +q 2 )/2m B . Further, one chooses the polarization vector of the D* meson 

as 

e(0) = — (bD.|,0, 0,E D .), e(±) = T^=(0,l,±t,0). (A.2) 

In this frame, we choose the polarization vector of the virtual gauge boson e, 
which can be, longitudinal (m — 0), transverse (m = ±) or timelike (m = t): 

e(0) = — j=(|p D (*)|, 0, 0, -g ) , e(±) = ^=(0, ±1, -i, 0) , 
V q V 2 

Kt) = ^= = -^(qo,0,0,-\p D *\), (A.3) 

The leptonic tensor will be evaluated in the q 2 rest frame. In this frame, we 
choose the transverse components of helicity basis e to remain the same and other 
two components are taken as 

e(0) = (0,0,0,-1), e(t) = (1,0,0,0). (A.4) 
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Let 9i be the angle between the three-momenta of D* meson and the charged 
lepton in the q 2 rest frame, and x be the opening angle between two decay planes. 
We define the momenta of the lepton and anti-neutrino pair as 



pf = (Ei, pi sin 61 cos x, pi sin 9[ sin x, —pi cos 9j) , 
Vu = (jPh ~pl sm Qi cos X) — pl srn @i sm X) pl cos @i) 



(A.5) 



where the lepton energy E\ = (q 2 + m 2 )/2\J~(f and the magnitude of its three- 
momenta is p\ = (q 2 — m^)/2yg^. 

A. 2 Form Factors 

The relevant form factors for the B — > D* matrix elements of the vector V„ = cry^b 



and axial-vector = 07^756 currents are defined as EC 



(D*{p D *,e*)\V,\B{p B )) = J %V }f} W^'PS. 



m B + m D * 



(D*(pD*) \B{jp B )) = 2m D *A {q 2 )—^ qil + (m B + m^A^q 2 ) 



C.q 

q 2 



-A 2 (q 2 



C.q 



(Pb + Pd* 



m 



B 



ni 



D* 



(A.6) 



In addition, from Eq. ( [A.6 ) one can show that the B — > D* matrix element for 
the scalar current vanishes and for the pseudoscalar current reduces to 



(D*(pn.,e*)\cr, B b\B(p B ) 



2m D *A (q 2 ) 
m b (fj) + m c (/i) 



e .q. 



(A.7) 



B. Form factors in the Heavy Quark Effective Theory 

In the heavy quark limit for the b, c quarks (m& )C ^> Aqcd), both charm and the 
bottom quark in the hadronic current have to be replaced by static quarks h v ' jC and 
h Vjb , where v B = p B /m B and v'p, = PD*/ m D* are the four-velocity of the B and 
D* mesons, respectively. The b — > c transition can be studied in the heavy quark 
effective theory (HQET). In this effective theory, the matrix elements of the vector 
and axial vector currents, and and A^ , between bottom and charm mesons [?] 
are defined as 

(D{v')\ \B{v)) = ^hl^B [h+(w) (v + v% + h^w) (v - v'), 
(D*(v f , e') \V lt \B(v)) = i^m B m D « h v (w) e'* u v' a , 



(D*(v',e')\ A„\B{v)) = ^m B m D « \h Al (w) (w + 1) e^* - h M {w) e 

-h A3 (w)e'*-vv' 



■vv„ 



(b.i; 
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where the kinematical variable w = vb-Vd* = (m B + m 2 Dt — q 2 )/2rriBmD* ■ 

Now we consider the form factors for the B — > D* matrix element in the HQET. 
The form factors h Ai (w) are related to the form factors in Eq. ( |A.6| ) fT3| , |32] , |35f in 
the following way, 

RoM 



T-, W + 1 1 t \ A I 2\ 

Rd* — —h Al (w), A [q ) 



R 



D* 



-h Al (w) , 



R2(W) 

Rd* 



h Al (w), V(q 2 ) 



Ri(w) 
Rd* 



h Al (w) 



(B.2) 



where R^* = 2^mBm* D / (itib + tti*d)- The w dependence of the form factors can be 
found in flTSL 



and the summary of the results are 



h Al (w) = h Al (l) 



1 - 8p 2 z + (53p 2 - 15)z 2 - (231p 2 - 9l)z 3 

R x {w) = Ri(l) - 0.12(w - 1) + 0.05(w - l) 2 , 
R 2 (w) = R 2 (l) + 0.11(w - 1) - 0.06(w - l) 2 , 
R ( w ) = R (\) - 0.11(w - 1) + 0.0% - l) 2 , 



(B.3) 



where z = (a/u> + 1 — \/2)/(y/w + 1 + \/2). The numerical values of the free param- 
eters p 2 , h Al (l), i?i (1) and i?2(l) are taken from p5|, 



^(1)1^1 = (34.6 ± 0.2 ± 1.0) x 10 
p 2 = 1.214 ±0.034 ±0.009, 
i?i(l) = 1.401 ±0.034 ±0.018, 
^(1) = 0.864 ± 0.024 ± 0.008 , 



-3 



(B.4) 



and Rq(1) = 1.14 is taken from Ul3 |. In the numerical analysis, we allow 10% 
uncertainties in the -Ro(l) value to account higher order corrections. 
In the HQET, the transversity amplitudes of Eq. ( [2.3|) become 



m B (l -n)(w + l)^/n , , 



l + 



{w-l)(l-R 2 (w)) 



r:-2nw) ' L (1-^*) 

A\\ =m B y/2r~*{w ± l)h Al {w){\ - g A ) , 
A± = - m B \/2r^(w 2 - l)h Al (w)Ri{w){\ + gy) , 



v 1 - 9a) , 



1 + r 2 — 2r*w 



■h Al (w)R (w) (1 - g A ) 



m 2 R (l ± r 2 - 2r* 



mi{m b {n) + m c {n)) 



9p 



(B.5) 



where r* = rriD*/mB, and 



A, 



(a± 



(B.6) 
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B.l Angular coefficients 

The expressions for the twelve angular coefficients in the B — > £)*(—>• Dir)l~isi 
angular distribution are summarized according to D* helicity combinations A1A2. 
The longitudinal V°'s (AiA 2 = 00) are given by 



^ = -2(l-^)|Xo| 2 , 



771 

V 3 ° = -S-LRe[A tP A* ] . 

The transverse V T, s (AiA 2 = ++, , H — , — h) are given by 

rn 1 ( m? s 



14 


r+ 


1^1 


\*\ 


i 2 + 


i^±i 



-4ReL4,i^ 



^ = -(l-^)(l^lll 2 -K| 2 ), 



V- 



r 



m. 



Im[A||A* 



The mixed \^ 0T 's (AiA 2 = 0±, ±0) are given by 

n 0T = V2(l-^)Re[An^], 

2 

V 2 0T = 2V2Re[- A ± A* + ^||A*p] , 
V° T = 2v^Im 

^ 4 0T = v^(l-^)lmL4 ± ^ 



(B.7) 



(B.8) 



(B.9) 
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